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Abstract 

We explain why aliasing can be detected in a generic temporally- 
sampled stationary signal process. We then define a concept of sta- 
tionarity that makes sense for single waveforms. (This is done with- 
out assuming that the waveform is a sample path of some underly- 
ing stochastic process.) We show how to use this concept to detect 
aliasing in sampled waveforms. The constraint that must be satis- 
fied to make detection of aliasing possible is shown to be fairly un- 
restrictive. We use simple harmonic signals to elucidate the method. 
We then demonstrate that the method works for continuous-spectrum 
signals — specifically, for time series from the Lorenz and Flossier sys- 
tems. Finally we explain how the method might permit the recovery of 
additional information about Fourier components outside the Nyquist 
band. 
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1 Introduction 



Detection of aliasing from temporal samples alone, with no restrictions on the 
original continuous-time source, is impossible because any set of samples may be 
reconstructed (using Shannon's sine filter) to a properly sampled signal having 
the same samples. However, quite often additional information about the source 
is available. It is, of course, obvious that tight constraints on the source would 
permit perfect reconstructions of vastly undersampled signals. For example, the 
constraint that the data comes from a linear function of time makes any two 
samples sufficient. A less extreme example is a signal with a lower as well as 
an upper frequency cutoff (a bandpass signal). For bandpass signals, it is well- 
known that one can sample at a rate below twice the highest frequency while still 
achieving perfect signal recovery (see Q p. 138, theorem 13.3]). 

What are the weakest constraints that one can put on the signal and still 
get something — detection of aliasing, for example? Here, we examine constraints 
of stationarity. In 1988 Hinich and Wolinsky || suggested a bispectral test for 
detecting aliasing in temporally sampled stationary stochastic processes^. The 
test aroused some controversy || Q which is examined in 0, 0. In JtJ we show, in 
detail, that the test does detect aliasing in some signal processes and that it is the 
constraint of stationarity that makes the detection of aliasing possible. Briefly, if 
we undersample a stationary process and then reconstruct a continuous-time signal 
from the samples using the Shannon sine filter, the reconstructed process will not, 
in general, be stationary. In contrast, a proper sampling followed by reconstruction 
will not destroy stationarity because this procedure just reconstructs the original 
signal. Detecting nonstationarity in the reconstructed process thus suffices to 
establish the existence of aliasing in the time series, provided it can be assumed 
that the original signal was stationary. These results are reviewed in Section |3[ 

Applying these concepts requires either a random sample of the paths of the 
process or an assumption of ergodicity which makes it possible to extract statistics 
from a single sample path. In this paper we attempt to generalize the results for 
stationary processes to the more common case where we have only a single sample 
path and can make no assumption of ergodicity. In other words, we look for ways 
to discover undersampling in a time series drawn from a single waveform, which 
may or may not be a sample path of some underlying stochastic process. We define 
sampling stationarity, a form of stationarity that makes sense for single waveforms, 
and show that it can be used to detect aliasing in complex, continuous-spectrum 
signals. We present reasons to believe that sampling stationarity should be a 

4 In the following, we will use the terms signal process or just process for stochastic 
signal processes. Except when we use the terms sample path for a realization of a stochastic 
process or random sample, the word "sample" will refer to temporal sampling. 
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generic^] property of signals and that the destruction of sampling stationarity by the 
process of undersampling and reconstruction should occur quite generally. Finally, 
we explain how it might be possible to use the reconstructed sample statistics plots 
(RSS plots) that we use to detect aliasing to obtain additional information about 
individual Fourier components beyond the Nyquist frequency. 

The remainder of this paper proceeds as follows. After illustrating the key idea 
of this paper with an example in Section^, we proceed, in Section ||], to demonstrate 
how a constraint of stationarity permits the detection of undersampling in some 
signal processes. Then in Section || we define sampling stationarity. In Section 41 



we use examples to show that the concept of sampling stationarity does, indeed, 
enable detection of aliasing for nontrivial signals. In Section we consider the 
case of periodic signals. For this class of signals, we provide a complete explanation 
of how (and when) the method of high-frequency detection works. The possible 
extension of this explanation to nonperiodic signals is then discussed in Section 4.3. 
In Section 4.4 , we present some reasons to believe that the plots that we have 
used to detect aliasing may also be used to recover some portion of the original 
signal's high-frequency content. This is followed by suggestions for further work 
(Section and a conclusion that summarizes the work in this paper (Section ||). 
Two appendices contain computational and mathematical details. 



2 Example 

The key idea of our approach is captured by a very simple example. Suppose that 
we sample a square wave that takes the values —1 and 1. There is a unique properly 
bandlimited signal that has this time series as its samples. We can compute this 
signal by applying the Shannon sine filter to our time series. We may regard this 
computation as an attempt to reconstruct the original continuous-time signal. If 
we can reject this reconstructed signal as the source of our samples, then we must 
conclude that the time series contains aliased components. 

Note that our given time series consists only of — l's and l's. The reconstructed 
signal, on the other hand, is necessarily a continuous function of time, taking on 
all values in the interval [—1, 1] (and, in fact, beyond). The only way that we could 
have obtained a sequence of —l's and l's by sampling the reconstructed signal is if 
we had chosen a particular sampling rate (or one of its subharmonics) and a unique 
shift of the sampling comb. Any other combination of sampling rate and shift 
would have produced a series that takes a continuum of values. The probability 

We use the term generic in a nontechnical sense. The term usually occurs in a situation 
where one would like to say "with probability 1" but where no obvious probability measure 
exists. 
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of having chosen the special sample rate and shift that give a sequence of — l's 
and l's is clearly zero, provided that our sampling rate was chosen independently 
of the source. With this proviso, then, we can reject (with probability 1) the 
hypothesis that our time series consisting of —l's and l's came from sampling the 
reconstructed signal. 

The assumption that the sampling rate was chosen independently of the source 
is justified in most (but not all) cases of practical importance because we can rule 
out any interdependence between the source and the sampling rate on physical 
grounds. For example, if a signal produced by a distant source is sampled at a 
predetermined rate, such a coupling is clearly out of the question — it would amount 
to believing that the process that produced the signal "knew" when we were going 
to sample at a distant location. 

We can draw valid conclusions from a sampled signal about Fourier compo- 
nents beyond the Nyquist frequency only if we can put constraints on the original 
continuous-time source. How can we characterize the constraints that we are im- 
posing in this case? Effectively, we are assuming that the sample times (which are 
determined by the sample rate and shift) do not play a distinguished role in the 
source. Showing that the sampling times are distinguished in the reconstructed 
signal then suffices to reject the reconstructed signal as the original source of the 
samples. 

How, then, can we extend this analysis to more general classes of signals? In 
the case of a square wave (or any signal that takes on a finite number of values), 
the appearance of the time series produced by sampling the reconstructed signal 
at the given sampling times could not be more different from the appearance of a 
time series produced by sampling at any other shift of the sampling comb. Thus it 
is clear what we mean when we say that the sample times are distinguished in the 
reconstructed signal. For more general signals, however, it is not so clear exactly 
what it means for the sample times to be distinguished. 

There is one obvious case in which we can be assured that the sample times are 
not distinguished in the original signal and in which we can detect the distinguished 
character of the sample times in the reconstructed signal. If the original signal is, 
in fact, a stationary signal process, then, by definition, no time is distinguished. 
The appearance of nonstationarity in the reconstructed signal would then indicate 
the presence of aliasing in the time series. The detection of aliasing in time series 
from stationary signal processes is the subject of the next section. Following that, 
we use our example of sampling from a square wave and insights from the case 
of stationary processes to develop a method for detection of aliasing in single 
waveforms. 
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3 Detection of Aliasing in Stationary Pro- 
cesses 



Consider the case of detection of aliasing in stationary signal processes. We start 
with the simplest stationary processes imaginable — randomly shifted periodic sig- 
nals. If we have a waveform, x(t), with period T, then we can produce a stationary 
process by adding to t a random time shift, 6, that is evenly distributed on [0, T). 
A sample path of our process then has the form x(t + 6) for a particular choice of 

e. 

Consider then the effect of undersampling and reconstruction on a simple sine 
process, 

x(t) = sin(27r/t + 27r/0), (1) 

where 6 is evenly distributed on [0,/ _1 ). If we undersample with a sampling 
interval At, corresponding to the Nyquist band [— (2Af) _1 , (2Ai) _1 ), and then 
reconstruct via convolution with the sine filter, we get the sine process given by 

X r (t) = sin(27T/t + 27T/0)- (2) 

Here, / is the aliased frequency, given by / = f + kf/At where kf is the unique 
integer that places / in the Nyquist band. The key point is that the phase of 
the reconstructed signal is the same as the phase of the source even though the 
frequency has changed to the aliased value /. For a process with a single harmonic, 
the reconstructed signal remains stationary because the phase term, 2nf8, is evenly 
distributed on 2ir. 

Consider then a second signal process, 

x(t) = sin(2vrat + 2irad) + sin(2vr/3t + 2ir/38), (3) 

where f3 is an integer multiple of a and 9 is chosen randomly from the interval 
[0,a -1 ). Since the time shift, 9, is the same for both components, this is, for the 
various values of 9, just a shifted waveform of a given shape. Since 6 is evenly 
distributed over the period, a -1 , the process is stationary. 

If we sample this process at a rate low enough for both components to be 
aliased and then reconstruct using the sine filter, we get 

x r (t) = sin(2vrat + 2-kolQ) + sin(2vr/3t + 2ir/36), (4) 

where a and (3 are the aliased frequencies. Although the phase terms, 2na9 and 
2irf38, are still evenly distributed over 2n, they now correspond to different time 
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Figure 1: Plot illustrating the nonstationarity of a randomly shifted waveform 
that has been undersampled and then reconstructed. The upper curve is the 
sixth moment^. The lower curve shows the process envelope. The original 
process is given by Equation |3] with a = 1.0, (3 = 3.0, 9 G [0, 1), and At = e. 

shifts for the two components. Thus, we no longer have a single shifted waveform 
and we can expect, in general, that stationarity will have been lost. 

We illustrate this loss of stationarity on an example by setting a = 1.0 and 
j3 = 3.0 in Equation || and choosing a sample time, At, equal to e. We may detect 
the loss of stationarity by examining the envelope of the reconstructed process. 
We define the envelope of a process, X t , as the support of the probability density 
of Xt as a function of t. For a process produced by randomly shifting a periodic 
waveform, the envelope may be conveniently displayed by plotting the sample 
paths corresponding to a representative collection of time shifts as in Figure |l|. 
Clearly, a stationary process must have a constant envelope. If we compute the 
envelope for the process defined in Equation |I| with the parameter values that 
we have specified, we get an oscillating figure (see Figure ||). This implies that 
the signal is nonstationary. In fact, it is cyclostationary with period equal to the 
sampling interval. 

As explained above, we do not lose stationarity when our original signal is a 
single sine wave. Nor do we lose stationarity when ratios between frequencies are 
preserved under the aliasing. For example, if At = 1.0 and the original frequencies 
are (10/9,20/9,30/9), they would alias to (1/9,2/9,3/9) and we would obtain 
another stationary process. But this situation is very special (non-generic). In 

7 The sixth moment was chosen for clarity of presentation. The second moment remains 
constant in this case. The fourth moment does oscillate but the scale of its oscillation is 
too small to allow meaningful display of the moment and the envelope on the same scale. 
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Figure 2: The sixth moment (upper curve) and the process envelope (lower 
curve) of the process given by Equation with a = 0.25, (3 = 0.75, 9 G [0, 4), 
and At = 1.0. 



general, more than one Fourier component is present and we do not have the 
special relationships between the sampling rate and the component frequencies that 
preserve ratios between frequencies when undersampling. Thus, we expect that, 
generically, a stationary process formed by randomly shifting a periodic waveform 
will lose stationarity upon undersampling and reconstruction. 

Within the context of single shifted waveforms, the destruction of stationar- 
ity can occur in some remarkable situations. Consider that undersampling and 
reconstruction can break stationarity even when only one of the components (/?, 
say) is aliased and when (5 aliases to a or —a. In other words, stationarity can 
be broken even when the two components, after undersampling, lie right on top 
of each other. We can see this by choosing a = 0.25, f3 = 0.75, and At = 1.0 
in Equation ||. The envelope for the reconstructed process is shown in Figure || 
where the nonstationarity is apparent. (Of course, we cannot possibly detect this 
loss of stationarity by examining only a single sample path, since the sample path 
will never be more than a single sine wave of some amplitude and phase.) 

Not all stationary processes are randomly shifted periodic waveforms. What 
can we say about more general stationary processes? It is clear that, if there exist 
no phase relationships between any of the Fourier components of the process, 
then undersampling and reconstruction will not destroy stationarity. For a generic 
stationary process, though, we would expect at least some sets of components to 
exhibit phase relations. In that case, we would expect stationarity to be destroyed 
because it is difficult to imagine how the destruction of stationarity associated with 
one set of components could somehow be canceled out by the presence of other 
incommensurate components. 
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This argument, together with the observation that there is simply no reason to 
believe that stationarity should be preserved under undersampling and reconstruc- 
tion, suggests that the loss of stationarity should be a general feature of stationary 
processes. 

4 Detection of Aliasing in Single Sample Paths 

The method that we have just used to detect aliasing in a sampled stationary 
process requires complete knowledge of the discrete-time process obtained by sam- 
pling the original continuous-time source. Usually, however, we have available to 
us only a single sample path. Therefore, we require a method for detecting aliasing 
in a single waveform which may or may not be a sample path of a stochastic signal 
process. 

We may develop such a method by reconsidering the example of sampling from 
a square wave discussed in Section || in light of our discussion of the effect of un- 
dersampling on stationary signal processes. Recall that the sampled time series 
from the square wave takes on the values —1 and 1. We may state this in statis- 
tical language by saying that the one-time probability density of the time series 
consists of two Dirac delta functions centered at —1 and 1, respectively. Now, we 
would have obtained the same one-time statistics if we had sampled the original 
square wave with any shift of the sampling comb. We will say that a waveform has 
sampling stationarity for a given sampling interval if the one-time sample statistics 
do not change as the position of the sample comb is shifted along the waveform. 
Observing that the original square wave had sampling stationarity for the given 
sampling interval is essentially equivalent to saying that the sample times were 
not distinguished in the source^ Note that the signal obtained by applying the 
Shannon sine filter to the time series does not have sampling stationarity — the 
one-time statistics of the reconstructed signal vary dramatically with shifts of the 
sampling comb (see Figure ^). This lack of sampling stationarity corresponds to 
the distinguished role of the sample times in the reconstructed signal. Of course, 
this distinguished role for the sample times is what allowed us to reject the recon- 
structed signal as a candidate for the original source of the samples and, thus, to 
conclude that the sampled series contained aliased components. 

8 Note that the original square wave does not have sampling stationarity for a sampling 
interval equal to its period. In general, a periodic signal will not have sampling stationarity 
with respect to sampling intervals commensurate with its period. However, the set of 
sampling intervals that are commensurate with a given period has Lebesgue measure 
zero. Clearly, the probability of choosing such a special sampling interval is under the 
assumption that the sample times are chosen independently of the source. 
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This discussion suggests the following test for aliasing in signals (no underlying 
stochastic process assumed). Collect the statistics on the recorded samples. Re- 
construct the signal at various shifts of the sampling comb and collect the statistics 
at these reconstructed samples. Compare with the original statistics. (We use the 
term "statistics" loosely, without the assumption that the samples are independent 
samples of some underlying probability distribution.) If we find that the recon- 
struction has different statistics at some shift of the sampling comb, an assumption 
of sampling stationarity for the original signal implies that the reconstruction is 
not the original signal and therefore that the signal was undersampled. 

For this test to be at all useful, two questions must be answered: 

1. Are typical signals characterized by sampling stationarity? 

2. Do typical undersamplings reconstruct to signals for which sampling sta- 
tionarity is violated? 

The answer to question p] is clearly "yes" for sample paths of ergodic stationary 
processes and for signals from ergodic dynamical systems. It is also clear that 
there are other classes of signals which possess sampling stationarity. For example, 
general periodic signals (not just square waves) possess sampling stationarity if the 
sampling interval is incommensurate with the signal period (a generic condition). 
Below, we conjecture that sampling stationarity is a generic property of signals. 

The examples that we consider next suggest that the answer to the second 
questions is also "yes" . 

4.1 Examples 

In each of the examples listed below, the time series to which we apply our test 
for aliasing was split into two interleaving series, D\ from the samples taken at 
[0,2Ai,4At, ...] and D2 from the samples taken at [Ai,3Ai,5Ai, ...]. The sample 
statistics corresponding to D\ and D2 are plotted in blue and green, respectively. 
For original signals with sampling stationarity, these two histograms will coincide. 
We then produce a reconstruction from D±, computed at the times corresponding 
to D2. The sample statistics corresponding to this reconstructed series are shown 
in red. If the red histogram is significantly different from the blue, then the 
reconstructed signal does not have sampling stationarity. 

Example: For a periodic signal, the generic condition of incommensurabil- 
ity of the sampling interval and the signal period implies that the signal has the 
property of sampling stationarity. But we also find that undersampling and re- 
constructing produces a signal that does NOT have sampling stationarity, as illus- 
trated in Figure ||. The data for the plot were generated by sampling (At = e/16) 
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Figure 3: Sample statistics of data and reconstruction from a square wave. 
The plot shows sample statistics of the data in blue and green (which are 
indistinguishable) and the reconstruction in red. (See the beginning of Sec- 
tion [O] for an explanation of the blue and green histograms.) The red 
histogram is obviously very different from the blue and green with which it 
would coincide if the reconstruction had sampling stationarity. The blue and 
green histograms have been rescaled so as to make the three histograms of 
comparable height. 



a sum of sines with frequencies (0, 1, 2, 10) and random amplitudes that ranged 
between .78 and 1.22. 

Example: If our signal consists of a sum of sine waves with incommensu- 
rate frequencies, then we cannot detect aliasing by this method (see Figure ||). 
Although such a signal will have sampling stationarity, the sampling stationarity 
will not be broken by undersampling and reconstruction because it is impossi- 
ble to have relationships between the phases of different Fourier components (see 
Section £D|). 

Example: The previous example might lead to the suspicion that this method 
works only for periodic signals (or step signals such as the square wave). However, 
the presence of incommensurate Fourier components does not necessarily destroy 
the ability to detect aliasing in a periodic waveform with more than one harmonic 
component. Figure || shows the result of combining a periodic waveform with 
incommensurate harmonics. The total power in the incommensurate harmonics is 
about 21% of the power in the periodic waveform. The sampling stationarity of 
the original signal and the breakdown of sampling stationarity with undersampling 
and reconstruction are apparent. This shows that, as long as some of our aliased 
Fourier components are commensurate with other components, the method can 
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Figure 4: Sample statistics of data and reconstruction from a periodic signal. 
The blue and green histograms coincide, indicating that the original signal 
had sampling stationarity. The red histogram, showing the statistics of the 
reconstructed signal, is obviously very different from the blue, showing that 
the reconstruction does not have sampling stationarity. 
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Figure 5: Sample statistics of data and reconstruction from a sum of sine 
waves with incommensurate frequencies. The blue and green histograms 
coincide. The red, showing the statistics of the reconstructed signal, is NOT 
obviously different. 
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Figure 6: Sample statistics of data and reconstruction from a mixture of 
a periodic waveform and incommensurate harmonics. The blue and green 
histograms coincide. The red, showing the statistics of the reconstructed 
signal, is obviously different. 

work. 

Example: So far, we have demonstrated that the method works for pure 
periodic signals and for periodic signals mixed with incommensurate harmonics. 
Figures [7|, ||, and ^ show that the method works for much more complex signals with 
continuous spectra. The signals are taken from the Lorenz and Rossler systems 
(see Appendix |SJ) . 

The success in detecting aliasing in time series from the Lorenz and Rossler 
systems suggests that the method may work for a very broad class of signals. 
Before attempting to determine how wide this class might actually be, we will 
look at the periodic case in order to begin to understand the precise mechanism 
of the method. 

4.2 A Closer Look at the Periodic Case 

If one samples a periodic signal incommensurately with the signal period, the 
samples end up mixing evenly around the waveform (see Appendix |B|). Thus, all 
shifts of a sampling comb with a sampling interval that is incommensurate with 
the period will produce the same statistics. This implies that: 

Theorem 1 A periodic signal will have sampling stationarity with respect to any 
sampling interval that is incommensurate with the period of the signal. 

Conversely, sampling with an interval that is commensurate with the period will, 
in general, produce statistics that depend on the sampling shift. Theorem [l] implies 
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Figure 7: Sample statistics of data and reconstruction from the x-coordinate 
of the Lorenz model. The blue and green histograms coincide and the red is 
clearly different. 
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Figure 8: Sample statistics of data and reconstruction from the x-coordinate 
of the Rossler model. The blue and green histograms coincide and the red is 
clearly different. 
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Figure 9: Sample statistics of data and reconstruction from the z-coordinate 
of the Rossler model. The figure speaks for itself. 



Theorem 2 Every periodic signal has sampling stationarity with respect to all 
sampling intervals except for a set of intervals with (Lebesgue) measure zero. 

Thus, the probability of choosing a sampling interval for which a given periodic 
signal does not have sampling stationarity is zero, provided the interval is chosen 
independently of the signal. 

What, then, is the effect of undersampling and reconstruction on this sampling 
stationarity? The sample statistics are determined by the shape of the waveform 
(see Appendix || for the exact formula). It can be shown that the process of 
undersampling and reconstruction is equivalent to sampling the original waveform 
at the same rate with the individual Fourier components shifted with respect to 
each other. When different components experience different time shifts, the shape 
of the effective waveform changes. Consequently, the statistics change. We now 
explain this in detail. 

When we sample a single harmonic with frequency / and phase (p every At 
time units, we get the values 



that: 



y n = sin (2irfnAt + ip) 



(5) 



We will temporarily suppress the phase and rewrite this expression as 



sin (2nfnAt) = sin (27r/nAt + 2-Kkn) 
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for any integer k, so that the reconstruction of this component at points 1 + s, 2 + 
s, 3 + s, ... is given by 



y n , s = sin 2tt / + -£ (n + s) At , (7) 



where the reconstruction chooses precisely one of the integral £;'s, which we will 
call kf, such that f + kf/ At is in the interval [—1/2 At, 1/2 At). We can now rewrite 
the reconstructed harmonic as 

Vn,s = sin (f + nAt + 2vr (f + |Q sAt 

= sin (2-KfnAt + 2itk f n + 2vr/sAt + 2vr£; / s) ' ^ ' 

= sin (2vr/nAt + 2vr/sAt + 27r/c/s) 

where we drop 2nkfn since kf is an integer. Thus, the reconstructed signal has 
samples at a shift, s, as though we were sampling the original waveform, but with 
the phase of the individual Fourier component shifted by the amount 2~KfsAt + 
2irkfS. The first term amounts to a time shift which is the same for all the 
components in the waveform. This implies that these first terms do not change 
the shape of the waveform and can be ignored. So we may consider the effective 
waveform (at a shift s) to be 

Aj sin(27r/jnAt + 2irk ft s + ifj) (9) 

i 

where we have reinserted the phase. The term 2irkf i s amounts to a time shift that 
is different for different /j. This difference in time shifts leads to a change in the 
shape of the effective waveform as s changes which in turn changes the sample 
statistics. 

For a given waveform, it is clear that almost any change in the shape of the 
waveform will change the sample statistics. (For example, a generic choice of s will 
change the heights of the extrema, changing the locations of the singularities in the 
histogram. See Appendix Jb[) Thus, we conclude that generically, periodic signals 
have sampling stationarity which is destroyed by undersampling and reconstruction. 



4.3 Nonperiodic Signals 

Now, we want to use the insight that we have gained for the case of periodic signals 
to get a better understanding of the answers to the two questions at the end of 
Section ||. We begin with some general questions about the kinds of signals to 
which our method might possibly apply. 
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Consider first the case of transient signals. In order to be able to talk about 
sampling stationarity at all, we have to be able to take as many samples as we 
want (at the given sampling rate) in order to be able to estimate the one-time 
probability distribution to arbitrary accuracy. This implies that we must think 
of our signals as functions of infinite time. In this context, any transient signal 
has trivial sampling stationarity — the probability distribution is a delta function 
at zero. By the same token, undersampling and reconstruction will not destroy 
this sampling stationarity. Thus, we need to restrict our attention to persistent 
(nontransient) signals. 

Within the class of persistent signals, it is clear that we need the signals that 
we consider to have well-defined sample statistics for arbitrary sampling intervals 
and shifts. Given that we are discussing aliasing, our signals also need to have 
a Fourier transform (in some sense). The set of signals with well-defined power 
spectra (which will have, in general, singular components) will clearly meet these 
criteria, although the actual class to which our method applies may be larger. In 
the following, then, we may take the term persistent signal to refer to a signal with 
a well-defined, nonzero power spectrum. 

Consider then the question of which signals have sampling stationarity for 
which sampling intervals. In the case of periodic signals, Theorems [l] and ^ provide 
what is essentially a complete answer — sampling stationarity holds for generic 
choices of signals and sampling intervals. At first glance one might try to generalize 
Theorem [l] to the following: 

Conjecture 1 Every signal has the property of sampling stationarity for every 
sampling interval At that is not commensurate with the period of any singular 
component of its spectrum. (FALSE) 

Unfortunately this conjecture is false as may be seen from the following coun- 
terexample. If we undersample and reconstruct a signal with a purely continuous 
spectrum (such as our signal from the Lorenz system), we will introduce no new 
singular components. Thus, the reconstructed signal will have a purely continuous 
spectrum. If the conjecture were true, then, such a reconstructed signal would 
have sampling stationarity for all sampling intervals by virtue of having no sin- 
gular components. Yet it is just the lack of sampling stationarity of this signal 
with respect to the given sampling interval that allows us to detect aliasing in this 
case. Thus, we know that there exist signals that lack sampling stationarity with 
respect to sampling intervals that are not commensurate with any singular compo- 
nent of their spectra and the conjecture is false. However, the reconstruction of an 
undersampled signal has a very special relationship to the interval with which the 
sampling was done. Thus, one expects that resampling the reconstruction with a 
new sampling interval not related to the original interval will yield statistics that 
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are again stationary with respect to shifts in the sampling comb. Therefore, we 
arrive at the following conjecture: 

Conjecture 2 Every signal has the property of sampling stationarity for every 
sampling interval At, except a set of At's with (Lebesgue) measure zero. 

This conjecture implies that, if one were to observe the reconstructed statistics 
varying with changes in the shift, this observation would be enough to conclude 
(with probability 1) that the samples came from an undersampled waveform. In 
other words, the truth of the conjecture would imply that the detection of under- 
sampling by the proposed method is generically free of false positives 

Next, we want to know when undersampling and reconstruction of persistent 
nonperiodic signals will yield new signals which have lost the property of sam- 
pling stationarity. (In other words, we also want to know when we can get false 
negatives.) The analysis that we have presented for periodic signals suggests that 
undersampling and reconstruction should destroy sampling stationarity for gen- 
eral persistent signals in which at least some of the aliased Fourier components are 
commensurate with other components of the signal. The reasoning is that each 
individual component can be regarded as a part of a family of harmonics and that 
the effective shape of the waveform associated with this family is changing with 
shifts of the sampling comb. There does not appear to be any reason to believe 
that combining different periodic waveforms, each of which is changing its sample 
statistics with shifts of the sampling comb, would result in sample statistics that 
do not change. Note that the condition that the original signal must have at least 
some commensurability will be satisfied by any signal with a nonzero continuous 
part to its spectrum as well as by periodic signals. Therefore, it seems likely that, 
generically, persistent signals have sampling stationarity that is destroyed by un- 
dersampling and reconstruction. In order to turn this last statement into a well 
defined conjecture, it will be necessary to define precisely what is meant by "gener- 
ically" in the case of persistent nonperiodic signals. The question of exactly how 
to define "persistent" must also be answered. Since the transformation that takes 
us from a waveform to sample statistics is extremely nonlinear, a proof is likely to 
be difficult. 

4.4 Recovery of High-Frequency Information 

The next question that presents itself is whether or not we can recover information 
about individual aliased Fourier components using the sampling-shift dependence 
of the reconstructed statistics. Ideally, we would like to know how much of the 
signal at an individual frequency, /, in the Nyquist band comes from each frequency 
that aliases to /. 
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Consider the one-time probability density of the reconstructed signal, p(x), 
as a function of both x and shift s. We call this two-dimensional surface a Re- 
constructed Sample Statistics (RSS) plot. The RSS plot has dependencies on s 
tied directly to the quantities kf i . Each in turn, determines the particular 
copy of the Nyquist band in which its corresponding is located. This chain 
of dependencies suggests that the RSS plot contains the information necessary to 
determine the contribution of each band to the signal at a given frequency in the 
Nyquist band. The inverse problem is greatly complicated by the interaction of 
the Fourier components and the nonlinear "projection" that turns the waveform 
into statistics. This extremely nonlinear inverse problem is the main object of our 
current research. 

5 Directions for Further Investigations 

In addition to the work already alluded to on the inverse problem formed by the 
RSS plots, there are other issues to explore. Included among them are: 

• What are the effects of noise on this method for detection of aliasing? 

• What is the effect of near commensurability of sample interval and signal 
period? 

• What is the effect of finite time-series length? 

• How does the departure of the statistics of the reconstructed signal from 
stationarity depend on the fraction of the total power that lies outside the 
Nyquist band? 

These questions are important to the practical usefulness of the method of 
high-frequency detection / recovery. 

6 Conclusion 

Although the idea of detection of aliasing is typically dismissed with references to 
the Nyquist criterion and the Shannon reconstruction theorem, we have demon- 
strated that detection of aliasing is possible with what appear, at first glance, to 
be very weak prior assumptions. The key concept is that of sampling stationarity. 
We emphasize that this concept makes sense for single waveforms. Although this 
concept arose in the consideration of step signals like square waves, its usefulness 
extends far beyond these signals. In particular, our method enables the detec- 
tion of aliasing in samples from nontrivial waveforms such as measurements from 
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motion on the Lorenz or Rossler attractors. As indicated above, many questions 
remain. Some of these are important for the practical utility of the concept of 
sampling stationarity and the associated RSS plots. 

A Computational Details 

The calculations represented in the paper were done with Matlab. The Lorenz 
equations, 

x =cr(y — x) 

y =x(R -z)-y (10) 
z =xy — bz, 

were integrated with parameter values of a = 10, R = 28, and 6 = 8/3 using 
Matlab's "ODE45" which is an adaptive step size routine. Relative tolerance was 
set to the default value of 1.0 x 10~ 3 and absolute tolerance was set to the default 
1.0 x 10~ 6 . Initial conditions were set at x = y = z = 1. Values for the x, y, and z 
coordinates were saved every 0.5 time units. 200,001 samples were taken and split 
into two interleaving time series each 100,000 samples long. The first series was 
used to reconstruct a signal via convolution with a sine filter of length 200,001. 
These very long series and filters were used to minimize the effects of truncating 
the convolution at the ends of the series. The histograms for the reconstructed 
signal were computed from the middle 50% of the reconstructed series. 
The Rossler equations, 

x = — z — y 

y=x + ay (11) 

z =b + z(x — c), 

were integrated in the same way with a sampling interval of 10 time units, and 
parameter values of a = 0.15, b = 0.2, and c = 10. In this way 200,001 samples 
were obtained and the splitting and reconstruction were done as described for the 
Lorenz equations. 

All the histograms presented here were originally calculated from much shorter 
time series (10,000 samples). The features that allow us to conclude that aliased 
components are present were all clearly visible in the histograms made from shorter 
time series although, of course, the histograms were considerably rougher. We 
conclude that the results that we have presented are certainly not an artifact of 
finite-length time series. 
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B Sample Statistics for a Periodic Signal 



If one samples a periodic signal, h(t), incommensurately with the signal period 
T, the samples end up mixing evenly around the waveform^]. The resulting his- 
togram is proportional to the reciprocal of the derivative of the waveform. This 
follows from the fact that the probability of getting any particular t (position along 
the waveform) is uniformly distributed over [0, T) which in turn implies that the 
probability of the interval [y, y + dy) is the probability of the corresponding dt or 
(l/T)(dy/h (t)). More precisely, the probability density for y is 

p(v) = y E (fc'&O)- 1 (12) 

t a eT(y) 

where 

T(y) = {t\te[0,T),h(t)=y} (13) 

Note that the density will have l/y/(y) singularities at the local maxima and 
minima of h(t). The form of the singularities follows from the fact that a generic 
waveform has maxima and minima with nonzero second derivative. 
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